Epidemic spreading on heterogeneous networks with identical infectivity 
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In this paper, we propose a modified susceptible-infected-recovered (SIR) model, in which each 
node is assigned with an identical capability of active contacts, A, at each time step. In contrast 
to the previous studies, we find that on scale-free networks, the density of the recovered individuals 
in the present model shows a threshold behavior. We obtain the analytical results using the mean- 
field theory and find that the threshold value equals 1/A, indicating that the threshold value is 
independent of the topology of the underlying network. The simulations agree well with the analytic 
results. Furthermore, we study the time behavior of the epidemic propagation and find a hierarchical 
dynamics with three plateaus. Once the highly connected hubs are reached, the infection pervades 
almost the whole network in a progressive cascade across smaller degree classes. Then, after the 
previously infected hubs are recovered, the disease can only propagate to the class of smallest degree 
till the infected individuals are all recovered. The present results could be of practical importance 
in the setup of dynamic control strategies. 

PACS numbers: 89.75.Hc, 87.23.Ge, 87.19.Xx, 05.45.Xt 



I. INTRODUCTION 



Many real-world systems can be described by com- 
plex networks, ranging from nature to society. Re- 
cently, power-law degree distributions have been ob- 
served in various networks Q, ■ One of the original, 
and still primary reasons for studying networks is to un- 
derstand the mechanisms by which diseases and other 
things, such as information and rumors spread over 0,0. 
For instance, the study of networks of sexual contact 
d H is helpful for us to understand and perhaps 
control the spread of sexually transmitted diseases. The 
susceptible-infected-susceptible (SIS) 0,0, susceptible- 
infected-removed (SIR) [Rj,[lll, and susceptible- infected 
(SI) |l2|, [l3|, [l4j models on complex networks have been 
extensively studied recently. In this paper, we mainly 
concentrate on the behaviors of SIR model. 

The standard SIR style contains some unexpected as- 
sumptions while being introduced to the scale-free (SF) 
networks directly, that is, each node's potential infection- 
activity (infectivity), measured by its possibly maximal 
contribution to the propagation process within one time 
step, is strictly equal to its degree. As a result, in the 
SF network the nodes with large degree, called hubs, will 
take the greater possession of the infectivity. This as- 
sumption cannot represent all the cases in the real world, 
owing to that the hub nodes may be only able to contact 
limited population at one period of time despite their 
wide acquaintance. The first striking example is that, in 
many existing peer-to-peer distributed systems, although 
their long-term communicating connectivity shows the 
scale-free characteristic, all peers have identical capa- 
bilities and responsibilities to communication at a short 



term, such as the Gnutella networks |15lll6| . Second, in 
the epidemic contact networks, the hub node has many 
acquaintances; however, he/she could not contact all 
his/her acquaintances within one time step |l7j . Third, 
in some email service systems, such as the Gmail system 
schemed out by Google , their clients are assigned by lim- 
ited capability to invite others to become Google-users 
[Tfij . The last, in network marketing processes, the refer- 
ral of a product to potential consumers costs money and 
time(e.g. a salesman has to make phone calls to persuade 
his social surrounding to buy the product). Therefore, 
generally speaking, the salesman will not make referrals 
to all his acquaintances pj| . Similar phenomena are com- 
mon in our daily lives. Consequently, different styles of 
the practical lives are thirst for research on it and that 
may delight us something interesting which could offer 
the direction in the real lives. 



II. THE MODEL 

First of all, we briefly review the standard SIR model. 
At each time step, each node adopts one of three possi- 
ble states and during one time step, the susceptible (S) 
node which is connected to the infected (I) one will be 
infected with a rate (3. Meanwhile, the infected nodes 
will be recovered (R) with a rate 7, defining the effec- 
tive spreading rate A = fi/j. Without losing generality, 
we can set 7=1. Accordingly, one can easily obtain 
the probability that a susceptible individual x will be 
infected at time step t to be 



X x (t) = l-(1-X) e ^ t - 1 \ 



(1) 



where 9{x, t — 1) denotes the number of contacts between 
x and its infected neighbors at time t — 1. For small A, 
one has 
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\ x (t)tt\6(x,t-l). 
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FIG. 1: (color online) R(oo) as a function of the effective 
spreading rate A on BA networks with (k) = 12, TV = 2000. 
The black line represents the case of standard SIR model, and 
the blue, green and red curves represent the present model 
with A = 4, 5 and 6, respectively. The arrows point at the 
critical points obtained from simulations. One can see clear 
from the inset that the analytic results agree well with the 
simulations. 



In the standard SIR network model, each individual will 
contact all its neighbors once at each time step, and 
therefore the infectivity of each node is equal to its de- 
gree and 9(x,t) is equal to the number of x's infected 
neighbors at time t. However, in the present model, we 
assume every individual has the same infectivity A. That 
is to say, at each time step, each infected individual will 
generate A contacts where A is a constant. Multiple 
contacts to one neighbor are allowed, and contacts not 
between susceptible and infected ones, although without 
any effect on the epidemic dynamics, are also counted 
just like the standard SIR model. The dynamical pro- 
cess starts with randomly selecting one infected node. 
During the first stage of the evolution, the number of in- 
fected nodes increases. Since this also implies a growth 
of the recovered population, the ineffective contacts be- 
come more frequent. After a while, in consequence, the 
infected population begins to decline. Eventually, it van- 
ishes and the evolution stops. Without special statement, 
all the following simulation results are obtained by aver- 
aging over 100 independent runs for each of 300 different 
realizations, based on the Barabasi-Albert (BA) [lj net- 
work model. 



III. SIMULATION AND RESULTS 

Toward the standard SIR model, Moreno et al. ob- 
tained the analytical value of threshold (k)/(k 2 ) [Tll |. 
Similarly, we consider the time evolution of I k (t), Sk(t) 



and Rk{t), which are the densities of infected, susceptible, 
and recovered nodes of degree k at time t, respectively. 
Clearly, these variables obey the normalization condition: 



h(t) + S k (t) + R k (t) = 1. 



(3) 



Global quantities such as the epidemic prevalence are 
therefore expressed by the average over the various con- 
nectivity classes; i.e., R(t) — ^ k P{k)R k {t). Using the 
mean-field approximation, the rate equations for the par- 
tial densities in a network characterized by the degree 
distribution P{k) can be written as: 



dh(t) 
dt 

dS k (t) 
dt 

dRkjt) 
dt 



.,,(, lt)ti ,(g gag ,, )) 



hit), 



(6) 



where the conditional probability P(k'\k) denotes the de- 
gree correlations that a vertex of degree k is connected 
to a vertex of degree k' . Considering the uncorrelated 
network, P(k'\k) = k'P(k')/{k), thus Eq. (4) takes the 
form: 

dIk[t) -hit) + ^ } S k (t)J2AP(k)Ik(t). (7) 



dt 



The equations (4-6), combined with the initial conditions 
R k {t) = 0, 4(0) = II, and Sfc(0) = completely de- 

fine the SIR model on any complex network with degree 
distribution P(k). We will consider in particular the case 
of a homogeneous initial distribution of infected nodes, 
I k = 1°. In this case, in the limit 1° —> 0, we can substi- 
tute I k (0) ~ and S k {0) — 1. Under this approximation 
and by taking the similar converting like from Eq. (4) to 
Eq. (7), Eq. (5) can be directly integrated, yielding 



S k (t) 



where the auxiliary function <j){t) is defined as: 



4>{t) 



l i: k AP(k)I k (t) _ ^ k AP(k)R k (t) 



io (k) (k) 

Focusing on the time evolution (j)(t), one has 

m) = j: k AP(k)i k (t) 

dt (k) 

J2 k AP(k)(l-R k (t)-S k (t)) 
(k) 

A am T. k AP{k)S k {t) 

- (k) m (k) 

A Y. k AP{k)e- Xk ^ 
(k) m (k) 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 
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Since ifc(oo) = and consequently limt—^ dcj)(t)/dt = 
0, we obtain from Eq. (13) the following self-consistent 
equation for c^oo as 



A 



— Xktfio. 



)• 



(14) 



The value (j)oo = is always a solution. In order to have 
a non-zero solution, the condition 



(k) dfioa 

must be fulfilled, which leads to 



A\ 



I 0oo=O 



=n> 1 



— J2kP(k) = XA>l. 



(15) 



(16) 
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This inequality defines the epidemic threshold 



1 

X 



(17) 



below which the epidemic prevalence is null, and above 
which it attains a finite value. Correspondingly, the pre- 
vious works about epidemic spreading in SF networks 
present us with completely new epidemic propagation 
scenarios that a highly heterogeneous structure will lead 
to the absence of any epidemic threshold. While, now, it 
is 1/A instead (see the simulation and analytic results in 
Fig. 1). Furthermore, we can also find that the larger of 
identical infectivity A, the higher of density of R(oo) for 
the same A from Fig. 1. 

From the analytical result, A c = 1/A, one can see that 
the threshold value is independent of the topology if the 
underlying network is valid for the mean-field approx- 
imation [20j. To further demonstrate this proposition, 
we next compare the simulation results on different net- 
works. From Fig. 2, one can find that the threshold 
values of random networks, BA networks with different 
average degrees, and BA networks with different sizes 
are the same, which strongly support the above analysis. 
Note that, in the standard SIR model, there exists obvi- 
ously finite-size effect while in the present model, 
there is no observed finite-size effect (see Fig. 3(c)). 



IV. VELOCITY AND HIERARCHICAL SPREAD 

For further understanding the spreading dynamics of 
the present model, we study the time behavior of the 
epidemic propagation. Originated from the Eq. (8), 
Sk(t) = e~ xk ^ l \ which result is valid for any value of 
the degree of k and the function <fr(t) is positive and 
monotonically increasing. This last fact implies that Sk 
is decreasing monotonically towards zero as time goes 
on. For any two values k > k , and whatever the ini- 
tial conditions S® and S®, are, there exists a time t' 
after which Sk{t) < Sk'(t). A more precise characteri- 
zation of the epidemic diffusion through the network can 
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FIG. 2: (color online) R(oo) as a function of the effective 
spreading rate A on BA and random networks with (k) — 6 
(a), BA networks for different attachment number m (m = 
(k)/2) (b), and BA networks with different size N (c). In 
figure (a) and (b), the network size is N = 2000, and in all 
the above three plots, the infectivity is A = 2. 
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that all of the curves show an initial plateau (see also a 
few previous works on hierarchical dynamics of the epi- 
demic spreading 0, l^, ), which denotes that the 
infection takes control of the large degree vertices firstly. 
Once the highly connected hubs are reached, the infec- 
tion pervades almost the whole network via a hierarchi- 
cal cascade across smaller degree classes. Thus, (ku(t)) 
decreases to a temporary plateau, which approximates 
(k) = 2m. At last, since the previously infected nodes 
recovered, all of which can be regarded as the barriers of 
spreading, the infection can only propagate to the small- 
est degree class. Then, the spreading process stops fleetly 
once the infected nodes are all recovered, as illustrated 
that {kn(t)) decreases to zero rapidly. 

Furthermore, we present the inverse participation ra- 
tio Yi(t) [24| to indicate the detailed information on the 
infection propagation. First we define the weight of re- 
covered individuals with degree k by Wk (i) = Rk (t) / ' R{t). 
The quantity Y 2 (t) is then defined as: 



Y 2 (t) = J2w 2 k (t). 



(21) 



Clearly, if Y 2 is small, the infected individuals are ho- 
mogeneously distributed among all degree classes; on the 
contrary, if Y 2 is relatively larger, then the infection is 
localized on some specific degree classes. As shown in 
Fig. 3(b), the function 5*2 (t) has a maximum at the early 
time stage, which implies that the infection is localized 
on the large degree k classes, as can also inferred from 
Fig. 3(a). Afterwards Y 2 (t) decreases, with the infection 
progressively invading the lower degree classes, and pro- 
viding a more homogeneous diffusion of infected vertices 
in the various degree classes. And then, Y 2 (t) increases 
gradually, which denotes the capillary invasion of the low- 
est degree classes. Finally, when Y 2 (t) slowly comes to 
the steady stage, the whole process ends. 



FIG. 3: (color online) Time behavior of the average degree of 
the newly recovered nodes (a) and inverse participation ratio 
Y 2 (b) in BA networks of size N = 5000. 



be achieved by studying some convenient quantities in 
numerical spreading experiments in BA networks. First, 
we measure the average degree of newly recovered nodes 
at time t, which is equal to the average degree of newly 
infected nodes at time t — 1, 



<M*)> 



J2 k k[R k (t)-Rk(t-l)} 
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1) 
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(18) 

(19) 
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In Fig. 3(a), we plot this quantity for BA networks as 
a function of the time t for different values of A and find 
a hierarchical dynamics with three plateaus. We can find 



V. CONCLUSION 

In this paper, we investigated the behaviors of SIR epi- 
demics with an identical infectivity A. In the standard 
SIR model, the capability of the infection totally relies 
on the node's degree, and therefore it leaves some prac- 
tical spreading behaviors alone, such as in the pear-to- 
pear, sexual contact, Gmail server system, and marketing 
networks. Accordingly, this work is of not only theo- 
retic interesting, but also practical value. We obtained 
the analytical result of the threshold value A c = 1/A, 
which agree well with the numerical simulation. In ad- 
dition, even though the activity of hub nodes are de- 
pressed in the present model, the hierarchical behavior 
of epidemic spreading is clearly observed, which is in ac- 
cordance with some real situations. For example, in the 
spreading of HIV in Africa , the high-risk population, 
such as druggers and homosexual men, are always firstly 
infected. And then, this disease diffuse to the general 
population. 



5 



Acknowledgments 

BHWang acknowledges the support of the National 
Basic Research Program of China (973 Program) under 
Grant No. 2006CB705500, the Special Research Founds 
for Theoretical Physics Frontier Problems under Grant 
No. A0524701, the Specialized Program under the Pres- 



idential Funds of the Chinese Academy of Science, and 
the National Natural Science Foundation of China under 
Grant Nos. 10472116, 10532060, and 10547004. TZhou 
acknowledges the support of the National Natural Sci- 
ence Foundation of China under Grant Nos. 70471033, 
70571074, and 70571075. 



[1] A. -L. Barabasi and R. Albert, Science 286, 509 (1999). 
[2] R. Albert, and A. -L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

[3] R. Pastor-Satorras, and A. Vespignani, Epidemics and 
immunization in scale-free networks. In: S. Bornholdt, 
and H. G. Schuster (eds.) Handbook of Graph and Net- 
works, Wiley- VCH, Berlin, 2003. 

[4] T. Zhou, Z. -Q. Fu, and B. -H. Wang, Prog. Nat. Sci. 16, 
452 (2006). 

[5] F. Liljeros, C. R. Rdling, L. A. N. Amaral, H. E. Stanley, 

and Y. Aberg, Nature 411, 907 (2001). 
[6] S. Gupta, R. M. Anderson, and R. M. May, AIDS 3, 

807-817 (1989). 
[7] M. MOrris, AIDS 97: Year in Review 11, 209-216 (1997). 
[8] R. Pastor-Satorras, and A. Vespignani, Phys. Rev. Lett. 

86, 3200 (2001). 
[9] R. Pastor-Satorras, and A. Vespignani, Phys. Rev. E 63, 

066117 (2001). 

[10] R. M. May, and A. L. Lloyd, Phys. Rev. E 64, 066112 
(2001). 

[11] Y. Moreno, R. Pastor-Satorras, and A. Vespignani, Eur. 

Phys. J. B 26, 521 (2002). 
[12] M. Barthelemy, A. Barrat, R. Pastor-Satorras, and A. 

Vespihnani, Phy. Rev. Lett. 92, 178701 (2004). 
[13] T. Zhou, G. Yan, and B. -H. Wang, Phys. Rev. E 71, 



046141 (2005). 
[14] A. Vazquez, Phys. Rev. Lett. 96, 038702 (2006). 



[15] M. A. Jovanovic, Modeling large-scale peer-to-peer net- 
works and a case study of Gnutell [M.S. Thesis], Univer- 
sity of Cincinnati, 2001. 



Http:/ /www.g 


;nutella.com/| 




T. Zhou, J. - 
Wang, arXiv: 


G. Liu, W. -J. Bai, G. Chen, and B. -H 
physics/0604083 (Phys. Rev. E In Press) 


Http:/ /mail.g 


Dogle.com /mail/help /intl /en/ about.html 



[19] B. J. Kim, T. Jun, J. Y. Kim, and M. Y. Choi, Physica 
A 360, 493 (2005). 



[20] Note that, if the connections of the underlying networks 
are localized (e.g. lattices), then the mean-field approxi- 
mation is incorrect and the threshold value is not equal 
to 1/A. 

[21] R. Pastor-Satorras, and A. Vespignani, Phys. Rev. E 65, 
035108 (2002). 

[22] M. Barthelemy, A. Barrat, R. Pastor-Satorras, and A. 

Vespignani, J. Theor. Biol. 235, 275 (2005). 
[23] G. Yan, T. Zhou, J. Wang, Z. -Q. Fu, and B. -H. Wang, 

Chin. Phys. Lett. 22, 510 (2005). 
[24] B. Derrida and H. Flyvbjerg, J. Phys. A 20, 5273 (1987). 
[251 W. -J- Bai, T . Zhou, and B. -H. Wang, arXiv: 

physics/0602173 



